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Quantum phase transitions in the one-dimensional period-two and uniform quantum compass 
model are studied by using the pseudo-spin transformation method and the trace map method. 
The exact solutions are presented, the fidelity, the nearest-neighbor pseudo-spin entanglement, spin 
and pseudo-spin correlation functions are then calculated. At the critical point, the fidelity and 
its susceptibility change substantially, the gap of pseudo-spin concurrence is observed, which scales 
as 1/JV (N is the system size). The spin correlation functions show smooth behavior around the 
critical point. In the period-two chain, the pseudo-spin correlation functions exhibit an oscillating 
behavior, which is absent in the uniform chain. The divergent correlation length at the critical point 
is demonstrated in the general trend for both cases. 
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I. INTRODUCTION 

Recently, the quantum compass model was introduced 
to describe some Mott insulators with orbit degeneracy 
by a pseudospinpl 0, where the coupling along one of 
bonds is an fsing type, but different spin components are 
active along other bond directions. The disorder effect in 
this model was also examined Q. The protected qubit is 
formed if it is separated from the low-energy excitations 
by a pseudo-spin excited gap. So a high quality factor, 
scalable and error-free scheme of quantum computation 
can be designed The symmetry of pseudo-spin Hamil- 
tonians is usually much lower than SU(2)[5[, and the 
result of numerical calculation has been shown that its 
eigenstates are at least twofold degenerate or highly de- 
generate and disordered @. The quantum XX-ZZ model, 
also called one-dimensional (ID) compass model, is con- 
structed by antiferromagnetic order of X and Z pseudo- 
spin components on odd and even bonds, respectively [7j. 
In addition, the ID quantum compass model is exactly 
the same as the ID reduced Kitaev model [9]. The ana- 
lytic eigenspectra in the latter model have been obtained, 
and it was shown that this model has one gapless phase. 
But the characters of the quantum phase transition have 
never been well studied previously. The realistic models 
of the orbital degeneracy are more complicated. 

For the compass model, the pseudo-spins may lead 
to enhanced quantum fluctuations near the quantum 
phase transitions (QPTs) and to entangled spin-orbital 
ground states. The numerical results have indicated 
that a first-order QPT occurs at J x = J z between 
two different states with spin ordering along either x 
or z directions 
delity 




, the ground-state (GS) fi- 
[l 4 \ \ M and entanglement 
emerged from quan- 
tum information science have been used in signaling the 
QPTs. To calculate these quantities accurately, it is nec- 
essary to know the exact GS wave function. The deriva- 



tives of the GS energy are intrinsically related to the GS 
fidelity l.'V. both can be used to identify the QPTs. For 
the special case of two spin —1/2 system, the entangle- 
ment is given by the concurrence. Quantum entangle- 
ment is one of the most striking consequences of quan- 
tum correlation in many-body systems, shows a deep re- 
lation with the QPT[l6|. Therefore understanding the 
entanglement is very important in QPTsJT?], EH- In the 
context of QPTs, the quantum entanglement have been 
the subject of considerable interests in the Dicke model 
[3 HE El and the XY model [HHl. 



On the other hand, experimental works on quasicrys- 
tals [25| and quasiperiodic superlattices |26| have in- 
spired theoretical interests in ID quasiperiodic systems. 
Period-two chain can be regarded as the intermediate 
one between uniform periodic chain and quasiperiodic 
chain, which have exhibited some unusual physical prop- 
erties. In this work, we study the one-dimensional com- 
pass model for both uniform and period-two cases by us- 
ing transfer matrix method[24j] and the method of Lieb, 
Schultz, and Mattis[27]]. The exact solutions for two cases 
are obtained. The GS fidelity and the energy gap be- 
tween uniform and period-two quantum spin chain are 
calculated. The behaviors of the pseudo-spin correlations 
with periodic boundary condition are given. 



The paper is organized as follows: In Section II, we 
give the model and the exact solution with periodic 
boundary condition. The calculation methods of fidelity 
and concurrence are introduced in Section III. The cor- 
relation functions are analyzed in Section IV. The paper 
is summarized in Section V, where we give some discus- 
sions and conclusions. 
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II. MODEL HAMILTONIAN AND EXACT 
SOLUTION 

The Hamiltonian of one-dimensional compass model is 
given by 

N' 

H = Y,iM°2i-i< + (1) 

i=l 

where J; is the nearest-neighbor interaction, are 
the Pauli matrix on site i, N = 2N 1 is the number of the 
sites, and (3 is the coupling parameter which determines 
the phase transition point. For Jn — J and J2i+i = a J , 
the model is a period-two case. By using the pseudo-spin 
(orbital) transformation method which is given by Brzez- 
icki et aZQ, we can define the modulated interactions for 
odd pairs of pseudo-spins {2i — 1, 2i} as — rf = erfi-i^f i' 
and the spin-flip operators of x direction are given by 
rf = ( — l)5Zfc=i s fc n^^ 2i crj. The two neighboring odd 
bonds can be expressed as the even {2i,2i + 1} bonds 
by a product —t[t[ +1 . Then the Hamiltonian of one- 
dimensional compass model can be written as follows 

N'-l 

H S = £[Ji(r?+/3T?7f +1 )] 

i=l 

-Jn>[t z n , + (2) 

Note that it looks like but is different from the transverse 
field Ising model. 

The vector s represents the state (si, • • • , sjv')- Here 
Si = 1 (si = 0) labels that the two pseudo-spins of 
the odd bond {2i — l,2i} are parallel (antiparallel). 

s = J2iLi s i 1S tne number of parallel odd pairs of spins. 
In this paper, we only discuss the ferromagnetic bound- 
ary condition of the quantum compass model, i.e. the 
case of the even s. The effective Hamiltonian (2) can 
be solved by using the Jordan- Wigner transformation for 
spin operators, 

7? = 1 - 2cJ Cl , (3) 
-rf = (ci + cj) n(l-2 C j Ci ), (4) 

j<i 

where Ci and cj are the anticommuting fermion operators. 
After this transformation, The effective Hamiltonian be- 
comes 

N'-l 

i=l 

Ji(3(ci4 +1 - cjc] +1 + c lCl+1 - cjc i+1 )] 
+J N '(3(cn'c\ - c f N ,c\ + c N >ci - c' N ,cx) 

N' 

+J N >c N ,c N i - Jj , (5) 



with 

ci = ci(-l) 1+s+ ^?=i c 3 c i, (6) 

Because we assume that the parity of s is even, it im- 
plies that only states with even numbers of Bogoliubov 
quasiparticles in the spectrum of the Hamiltonian (5). 
Under the periodic boundary condition (cat'+i = ci), the 
number of c fermions must be odd parity, as can easily 
be obtained from equation (6). Then the general form of 
the Hamiltonian is simplified to 



N' 

H g = J2i 2j i4 c i + Jip(ci4 

z=l 
N' 

t=i 



- 44+i + c * c *+i - 4 c *+i)] 

(7) 



For the period-two case, we can rewrite (7) as the follow- 
ing form by neglecting the last constant term, 

N> 

F = E + 2^ Bi ^ + h -°- )] ' (8) 

where the nonzero elements of the matrices A and B are 
given by 

Aij = 2JiSij — JifiSj^+i — Jj/3Sj t i—i, 

Bij — —Jif35j.i+i + JjPSjj-x; 
A\n' — An'i — —Jn<P, 
Bin' = —Bn'i = Jn'P- 

Equation (8) can be diagonalized by using the Bogoliubov 
transformation 

1 N ' 

Vk = r ^ii^kt + ^ J ki) c i + - 1pki)4)' 

i=l 

1 N ' 

% = g + ^i)4 + (<f>ki - ^ki)ci], (9) 

i=l 

where ipki is the eigenvector of the matrix (A+B)(A — B) 
and 4>ki is that of the matrix (A — B)(A + B). The 
eigenvalues of both matrices are corresponding to A^. We 
take k = 0, ±^,±2^7, This relation is satisfied 
with the periodic boundary condition. In general, the two 
eigenvectors {4>ki and tpki ) satisfy the following equations 

(A - B)$ k = A k $ k , (A + B)$ k = A fc &, (10) 

where <f> k and ip k are two column vectors. The diagonal- 
ized result takes the form 

H = J2 A Mvk-\). (11) 
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The excitation energies > 0. At zero temperature, the 
QPT points are those parameters that satisfy the condi- 
tion Afc = 0, and the two coupled coefficients of the Bo- 
goliubov transformation satisfy the following equations: 



Afe^fe, 
Afei/'fe, 



2JiP4>k,i 



(12) 



which can be derived from equation (10). For the period- 
two case, i.e. Jn — J and Ji%+\ = Ja, if we take 
J = 1 and assume that Vfe,2n = Ae l2nk and tpk,2n+i = 
g e i(2n+i)k ^ exact results of A^ can be obtained an- 
alytically from the coupled equations (12) by using the 
trace map method. The result is expressed as 



A 2 
A fe± 



= ±V4J 4 (/3 2 + l) 2 (a 2 - l) 2 
+ (a 2 + 1)(2J 2 /? 2 + 2J 2 ). 



64a 2 J 4 P 2 cos 2 fc 
(13) 



The excitation energies have two branches (A^_ and 
Afc-|_). For a special case a = 1, i.e., the uniform pe- 
riodic chain, the excitation energies can be simplified as 
2Jy/l + f3 2 — 2/?cosfc, which is the same as that in Ref. 
Q ■ The QPT point is determined by Aj,_ . At the critical 
point, the equation can be decoupled for A& = 0. One 
of equation (12) is rewritten as (j>k,i+i — \4>k,%- Due to 

the periodic boundary condition, (4) =1 should be 
satisfied. The only possibility is j3 = 1, i.e., there is only 
one QPT point at f3 = 1 in this case. On the other hand, 
the GS energy is expressed as E = — i ^ fe At which in- 
cludes the spectra of the ± branches. In thermodynamic 
limit, the summation can be replaced by an integral 



En 



-JN'±- ( (A fc _ +A fc+ 
27r Jo 



)dk. 



(14) 



The pseudo-spin excitation gap A, which is energy differ- 
ence between the first excited state and the ground state, 
is equal to Ao_ , which disappears at /3 = 1. 

From the Fig.l, we can find that the symmetries 
of the pseudo-spin gaps are broken more obviously 
as [3 is away from the QPT point in the period- 
two model. The symmetries remain for the uniform 
model. The quantum critical point is fixed at j3 c = 1 
which separates the disorder phase. In the vicinity of 
the quantum critical point, the linear relation A = 

y/ 10(a 2 + 1) - 2V25a 4 + 14a 2 + 25|1 - (3\ is generally 
satisfied. 



III. FIDELITY AND PSEUDO SPIN 
CONCURRENCE 

The exact GS wave function of the system must be 
obtained in order to calculate the fidelity and concur- 
rence. Similar to the Bardeen, Cooper, and Schrieffer GS 
wave function, we can write the present GS wavefunction 

as 



28] 



l*oC9)) =l[vk\Vac) for all k. 



(15) 



ax 




FIG. 1: Pseudo-spin excitation gap A on uniform and period- 
two cases of the compass model. The gaps collapse at the 
quantum phase transition point at j3 = 1 for different values 
of a. 



According to equation (9) and the definition of the 
fidelity [n( 



F(/9,<5) = I <*o(/3)|*oG0 + *)>|, 



(16) 



where 8 is a small quantity (5 = 10 is taken in our cal- 
culation), the fidelity and its susceptibility can be given 
by 



21im 1 -y>. 

s^o S 2 



(17) 
(18) 



The numerical results for the GS fidelity and its suscep- 
tibility are plotted in Fig. 2. An abrupt jump occurs in 
the vicinity of the QPT point (/3 C = 1) as a consequence 
of the dramatic change of the structure of the GS. It 
agrees exactly with our analytical derivations. One can 
see level-crossing at f3 = (3 C , indicating the first-order 
QPT in this model. 

In recent years, the concept of concurrence is usually 
adopted as the measure of the entanglement in spin —5 
systems. We will give the nearest-neighbor pseudo-spin 
two-point correlation functions to calculate the nearest- 
neighbor concurrence (NNC) of the system. Because of 
the reflection symmetry, the global phase flip symmetry, 
and the Hamiltonian being real, the nonzero elements are 
given by [H, H3] 



•i+ll 



'i+ll 



Gi,i+1, ( T i T i+l) — Gi+l,ii 

Gi t iGi + i t i + i — Gi y i + iGi + i t i, 
G4 _ j . 



(19) 



where Gj » = —^2 k il>ki4>kj- The definition of concurrence 
max[ri(i,j) - r 2 (i,j) - r 3 (i,j) - 
are the square roots of the 



is given by C{i,j) 
r±{i,j),0\, where 
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FIG. 2: The fidelity and the susceptibility of the period- two 
compass model versus (3 for a = 0.2,0.5,0.8,1.0 and N' — 
100. The first QPT point is obviously found at (3 C = 1. 




FIG. 3: The concurrence C versus (3 for a = 0.2,0.5,0.8, 1.0 
(N' — 100). The inset shows the derivative dpC as a function 
of (3. 



FIG. 4: The NNC as a function of (3 with a = 0.6 for N' = 
100. A gap AC for the concurrence is found at the critical 
point. The inset shows the size scaling of the gap. 



even-pseudo-spin chains with periodic boundary condi- 
tion in this paper. The QPT in the ID compass model 
is of first-orderpl 0] , the scaling behaviors at the critical 
point should be absent. But the discontinuousness of the 
concurrence at the QTP may exhibit the finite-size scal- 
ing behavior N — > — 1 [2^| , consistent with the present 
observation. Due to the concurrence gap, the value of 
dffC becomes minimum at the critical point. However, 
the maximum value of the concurrence occurs below (3 C 
is not related to the critical point. The present results 
for the concurrence are similar to those in the periodic 
quantum Ising chain model [l8|. 



eigenvalues of the product matrix R = pijfkj in de- 
scending order. The spin flipped matrix plj is denned 
as fa = (av <g> a v )p* J (a v <g> er a ). The p %j is the den- 
sity matrix for a pair of qubits from a multi-qubit state. 
In this way we can calculate the NNC of pseudo-spins. 
For the period-two chain, the concurrence C2i.2i+i and 
C2i+i,2t+2 are different. So we use the average concur- 
rence C = 5(C 2 i,2i+l + C 2 i+l.2i+2)- 

The numerical results for the concurrence as a function 
of (3 are given in Fig. 3. It is shown that the maximum 
value of the concurrence gradually increases with the in- 
crease of parameter a. If a is small enough, the entan- 
glement of nearest-neighbor pseudo-spins disappears in 
the larger j3 regime. A cusp of the first derivative of the 
concurrence occurs at the critical point (3=1, similar to 
those in Ref. [161 ]. 

A gap is found in the curve of NNC versus (3 at the 
QPT point (3 C in our calculation of the pseudo-spin con- 
currence. If the pseudo-spin chain goes to infinite, the 
gap has the critical behavior with AC oc N" 1 , as shown 
in the inset of Fig. 4. Obviously, it is the finite-size ef- 
fect. The question then arises: what is the origin of 
the concurrence gap? The answer is the symmetry of 
system which has been assumed by the ferromagnetic 



IV. SPIN AND PSEUDO-SPIN CORRELATION 
FUNCTIONS 

Firstly we show the numerical results of the ground- 
state spin correlations on odd {2i — 1, 2i} and even 
{2i, 2i + 1} bonds as a function of (3 with a periodic 
boundary condition. The value of < cli-i^fi > grad- 
ually increases with (3 while < crfjCrfj +1 > decreases 
with (3, as shown in Fig. 5. The crossing points of 

< c r | i _ 1 cr| i > and < cr^cr^i+i > curves for the same a 
occur at the quantum critical point. Actually, the com- 
pass model is a kind of pseudo-spin Ising chains at (3 = 
and (3 — > oo. As a result, the curves of spin correla- 
tions versus (3 are asymmetric. So < cr^i-i^ii >~ * and 

< crfj(T| l+1 >^ — 1 as (3 — > oo. 

It is found that the correlation gradually increases with 
the decreasing a. When a = 1, the numerical result at 
the critical point is the same as the analytical result by 
Brzezicki et al. 0j- 

Finally, we calculate the distance dependence of the 
pseudo-spin correlator < Tfrf +r > under the periodic 
boundary condition for the period-two and uniform cases. 
The the two-point correlation function is given bv[27| 
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FIG. 5: Spin correlation functions in the period-two chain 
and uniform chain for a = 0.2, 0.6, 1.0 and N = 100. 
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FIG. 6: Distance dependence of < TiT? +r > correlator at the 
critical point. The parameters are a = 0.6, 1.0 and N' = 200. 



to find that the oscillation occurs for a ^ 1, i.e. for 
period-two chain, which can be attributed to the different 
coupling coefficients of odd and even bonds. However, 
the similar trend appears in both cases, as shown in Fig. 
6. 



V. SUMMARY AND DISCUSSION 

By using the pseudo-spin transformation method and 
the trace map method, we obtain the exact solution of 
one-dimensional compass model with periodic boundary 
condition. The parameter a determines the symmetries 
of finite pseudo-spin excitation gap A, but the phase 
transition point is still fixed at j3 = 1. The quantum crit- 
ical point separates the disorder phase. The pseudo-spin 
liquid disordered ground state is the universal features 
in the ID compass model. The numerical methods to 
calculate the fidelity and concurrence are also given. We 
observe a first-order quantum phase transition between 
two different disordered phase. The concurrence gap AC 
displays the scaling property N = — 1. The spin and 
pseudo-spin correlation functions are calculated. Curves 
for the two spin correlation function cross exactly at the 
critical point for any value of a. It is observed that the 
distance dependence of < Tfrf +r > correlator displays 
oscillation in the period-two case, and a divergent cor- 
relation length at the critical point is observed in both 
uniform and period-two chains. 
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T i+r 



Gi+l,i+l 



G 



i+r-l,i+l 



G 



i+r-l,i+2 



G 



i+r- 



-l,i+r 

(20) 

which has the form of Toeplitz determinant. When 
r —>■ co, the correlators gradually decrease and approach 
the asymptotic value for large r in an algebraic wayQ. 
This correlator is positive for all r, indicating that there 
is the long-range ferromagnetic order. It is interesting 
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